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ABSTRACT 

Aims. We present a determination of the intrinsic colors and ages of galaxies at very high redshift, in particular old galaxies (OGs) 
within extremely red objects (EROs). To date, the definition of EROs has been restricted to objects with z < 2.5, however, here 
consider objects with the same properties but shifted to higher redshifts (z > 2.5). We therefore, refer to these objects as very high- 
redshift EROs (Z-EROS, herein). 

Methods. We analyze 63,550 galaxies selected in the XMM-LSS field. To obtain a reasonably sized sample of EROs, it is essential to 
consider a very wide area surveys. We identify targets within an area of 0.77 square degrees for which optical to mid-infrared data are 
available from SUBARU, UKIDSS, and Spitzer. We select Z-EROs based on their colors, and then perform a selection of only OGs. 
One of our novel innovations is to adapt the traditional method of EROs selection based on the filters I and K, to higher redshifts. 
Using our method, we identify 20 objects that satisfy the conditions required to be Z-EROs/OGs at redshifts 2.5 < z < 3.8 within 
some conservative constraints of errors in the photometric redshifts. For each of these objects, we calculate the corresponding color 
at rest (B-V), and estimate their average stellar mass and age by comparing this color with a synthesis model. 

Results. Our selected galaxies have high stellar masses (~ 10" M Q ) and are older than 1 Gyr, hence their stellar populations were 
formed at z > 4.7. After including additional galaxies with z < 2.5 analyzed in a previous paper, we find that the formation epoch 
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the stellar mass, the lower the age of the Universe at which it was formed. This result appears to conflict with A-CDM models that 
claim that the most massive galaxies formed after lower mass. 

Key words, galaxies: evolution — galaxies: high redshift — galaxies: statistics — galaxies: formation — Infrared: galaxies 



An important goal of extragalactic astronomy is to understand 
how stars are assembled into galaxies and how this assembly is 
related to their evolution. In hierarchical A-CDM models, star 
formation starts out in low-mass systems, which build more mas- 
sive galaxies through sequential merging (e.g., White & Rees 
1978; Somerville 2005; De Lucia et al. 2006). In this picture, 
the most massive galaxies should be found at relatively low red- 
shifts. This is apparently at odds with observations of massive 
galaxies at high redshift: a significant population of galaxies 
with stellar mass ~10 n M has been found at z~ 2- 3 (e.g., 
Fontana et al. 2004; Yan et al. 2004; Daddi et al. 2005; Rudnick 
et al. 2006; van Dokkum et al. 2006; Cassata et al. 2008; Kriek 
et al. 2008, 2009; Hempel et al. 2011; Lopez-Corredoira 2010 
[hereafter L10]). Stellar population synthesis models combined 
with broadband photometric data have been used to show that 
many of these galaxies contain an old stellar population, with 
ages that indicate a star formation phase occurred within 1-3 Gyr 
of the Big Bang. Even at redshifts higher than 3, there is evidence 
of massive early-type galaxies (Toft et al. 2005; Chen & Marzke 
2004; Rodighiero et al. 2007; Wiklind et al. 2008). These in- 
triguing results have motivated efforts to find out more about 
these high redshift massive galaxies, particularly the epoch at 
which they formed, which provides an important check of their 
compatibility with cosmological models. 

It is also often accepted that star formation in galaxies was 
more active at high redshift and therefore that these galaxies are 
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in general intrinsically bluer than at low redshift (e.g., Dickinson 
et al. 2003; Rudnick et al. 2006; Labbe et al. 2007), as would 
be expected if their populations were younger and had lower 
mass/luminosity ratios. However, there are also very red galax- 
ies at very high redshift, it is indeed well-established that at 
intermediate-high redshifts a significant population of red pas- 
sively evolving early-type old galaxies" (OGs) can be found in 
the field among extremely red objects (Cimatti et al. 2002; Yan et 
al. 2004; see McCarthy 2004 for a review), together with "dusty- 
reddened starbursts" (DSs). All these galaxies are "extremely red 
objects" (EROs). In our case, we are interested not in DSs but in 
finding old galaxies at high redshift. The analysis of EROs/OGs 
provides a convenient means of achieving this aim because we 
can derive their ages using simple synthesis models. Since EROs 
are a mixture of OG and DS populations, it is necessary to re- 
move DSs from ERO samples before the OG population can be 
investigated (Miyazaki et al. 2003, Fang et al. 2009). 

In this paper, we continue the studies initiated by L 10 to de- 
terminate the intrinsic color and age variations at different red- 
shifts in a statistical way for very red passively evolving galax- 
ies. While L10 analyzed EROs/OGs with z < 2.5, we analyze 
EROs/OGs at higher redshift (z > 2.5). Since the EROs are rare 
and clustered (Daddi et al. 2000; Roche et al. 2002; Miyazaki 
et al. 2003; Georgakakis et al. 2005), and even rarer at z > 2.5, 
wide-field surveys are essential to help us study their statisti- 
cal properties, hence, in this paper we use the large area survey 
XMM-LSS. The paper is organized as follows. In ^2) we de- 
scribe the criteria used to find the candidates to very red galax- 
ies. In ^3] we describe our sample and photometric data. In ^4] 
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we calculate the colors at rest of the galaxies. In §0 we derive 
the ages of these galaxies. And in §|6] we discuss the results. 
Throughout this paper, we assume a A-CDM cosmology with 
Q m = 0.24, Q. A = 0.76, and H = 73 Km/s/Mpc. 

2. Criteria for selecting extremely red objects/old 
galaxies at z > 2.5 

It is possible to use colors to identify EROs/OGs within a sample 
of galaxies. For instance, among the many references, Fang et al. 
(2009) and Pozzetti & Mannucci (2000) presented methods to 
select EROs/OGs in the UDF field using magnitudes in the three 
filters z'7750, J, and K. The criterion for selecting EROs is (using 
always AB magnitudes; Fang et al. 2009). 

(i 7150 - K) > 2A2, (1) 

and likewise, the criterion to select OGs within EROs is 

(J - K) < 0.20(17750 ~K) + 0.39, (2) 

which is valid only for galaxies with redshifts 0.8 < z < 2.5. 
There are also other methods for selecting OGs within ERO 
galaxies, and Fang et al. (2009) demonstrated that all of them 
are more or less consistent with each other. There may be some 
incorrect identification of OGs among the galaxies selected with 
this color method (Fang et al. 2009 estimated the contamination 
to be < 25%), and consequently may be present dusty galax- 
ies in the sample selected using this criterion, but the statistical 
comparison with the SED fitting method in Fig. 1 of L10 shows 
that there are neither significant differences nor systematic ef- 
fects with redshift. 

Equations ([1) and (f2]i for a galaxy at z — 1 .5 (the approximate 
center of the redshift interval) are equivalent to 

/ d log F v \ Id log F v \ 

\^) (J ' K)<0M+0 ^-d^) (h ^ (4) 

where the ((...))(filter 1, filter 2) expressions indicate the average 
over the wavelength range between the center of filter 1 and the 
center of filter 2. Moreover, if we translate these equations in to 
the equivalent one at higher redshifts (taking the closest AB filter 
redshifted instead of filters Z7750, J, K), we get 

J -[3.6] >2.55, 

K - [3.6] < 0.39 + 0.19(7 - [3.6]), 2.5 < z < 3.5 
H- [4.5] > 2.41, 

K - [4.5] < 0.56 + 0.29(7/ - [4.5]), 3.5 < z < 4.5 
H - [5.8] > 3.02, 

[3.6]-[5.8] < 0.35 + 0.14(7/ - [5.8]), 4.5 < z < 5.5 ' (5) 
#-[5.8] >2.43, 

[3.6]-[5.8] < 0.35 + 0.18(AT - [5.8]), 5.5 < z < 6.5 
#-[8.0] > 3.21, 

[3.6]-[8.0] < 0.59 + 0.23(K - [8.0]), 6.5 < z < 7.5 

where [x] stands for the magnitude at filter of x /mi. These are 
the conditions the EROs/OGs sources must satisfy when red- 
shifted to z > 2.5. We refer to these sources as Z-EROs/OGS. 



We apply this selection to more than 60,000 galaxies, as we see 
below. For a more secure selection and to avoid contamination 
by starbursts, we also apply another condition on the ratio of 
fluxes in K to [24.0] (Fontana et al. 2009) given by 

F([24])/F(K) < 6. (6) 

This criterion is valid for redshifts of 2.5-4.0. As explained in 
Fontana et al. (2009), for quiescent and "red and dead" galaxies, 
this criterion effectively applies a threshold to the specific star 
formation rate by comparing with theoretical predictions. 

3. Data 

Our main purpose is to extend the sample of L10 to include 
galaxies at z > 2.5, which were not analyzed by L10. We re- 
quire flux data for at least five filters from B to [8.0], including 
always the K filter and then three appropriate filters in the range 
required for the two color determination according to the sys- 
tem of equations given in Eq. ((5), with a signal-to-noise ratio 
(S/N) greater than 3. We also insist that there is a small amount 
of emission in the filter [24.0] to apply the selection criterion of 
Eq. ©. If we do not have any signal in that band, we assume 
that the object has no detectable emission at [24.0]. 

After exploring the data of the different available public sur- 
veys, we realized that we needed to combine a mid-infrared sur- 
vey with both optical and near-infrared ones covering large areas 
and deep enough to reach the high redshift sources. The survey 
UKIDSS (Hewett et al. 2006; Lawrence et al., 2007) is perhaps 
the highest quality survey in near infrared with these characteris- 
tics. It consists of five surveys of complementary combinations 
of depth and area, employing the wavelength range 0.83-2.37 
/mi in up to five filters ZY JHK, and extending over both high 
and low Galactic latitude regions of the sky. The UKIDSS uses 
the Wide Field Camera mounted on the UK Infrared Telescope 
(UKIRT) (Casali et al. 2007). Our study uses data from the DR3 
(Data Release 3) of the Ultra Deep Survey (UDS), specifically 
the XMM-LSS field centered on RA=02hl8m00s and DEC=- 
05d00m00s (J2000), which is a deep, wide survey with an area 
of 0.77 deg 2 and a 5<r depth up to J ~ 24.9, H ~ 24.2, and K ~ 
24.6 in the AB system (Lawrence et al. 2007). An early data re- 
lease of UKIDSS was previously used to study EROs (Simpson 
etal. 2006). 

The optical imaging observations of this field were car- 
ried out using the prime-focus camera (Suprime-Cam; Miyazaki 
et al. 2002) on the Subaru Telescope in the period from 
2002 September to 2005 September as part of the project 
S UB ARU/XMM-Newton Deep Survey (SXDSS; Furusawa et al. 
2008). The layout of the pointings is arranged in a cross shape so 
that each of the north-south and the east-west directions has an 
extent of ~1.3 deg (1.2 deg 2 ). This corresponds to a field span 
of a transverse dimension of ~75 Mpc at z~ 1 and -145 Mpc 
at z~ 3 in the comoving scale for (h,Q.M,Q\) =( 0.7,0.3,0.7 ). 
Observations of the SXDF were performed using five broadband 
filters, B, V, R c , i', and z' to cover the entire wavelength range 
observable with Suprime-Cam. The limiting magnitude in each 
band is B=28.2, V=27.8, R c =27.5, i'=27.2, and z'=26.3 (AB, 5 
cr, (p = 2" ) in the deepest images of the five pointings. 

We have the additional advantage that one of the regions 
covered by UKIDSS and SUBARU was covered by the Spitzer 
Wide-area InfraRed Extragalactic survey (SWIRE; Lonsdale et 
al. 2004). The SWIRE XMM-LSS field survey has an area of 9.1 
deg 2 in all seven Spitzer bands: Infrared Array Camera (IRAC; 
Fazio et al. 2004) at 3.6, 4.5, 5.6, and 8 /mi; and Multiband 
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Imaging Photometer (MIPS; Rieke et al. 2004) at 24, 70, and 
160 fim. The sensitivities in all the bands are 3.7 fiiy, 5.4 fiJy, 48 
jjjy, and 37.8 fi]y for the IRAC bands and 230 fiJy, 18 mJy, and 
15.5 mJy for the MIPS bands. The Spitzer Space Telescope has 
permitted the selection of objects in a spectral region (3-8 fim) 
that was not previously accessible from the ground. The selec- 
tion of objects at these wavelengths is particularly invaluable for 
high redshift massive galaxies because a 3-8 fim selection allows 
a sampling of the rest-frame near-IR flux of the high-z galaxy 
spectral energy distributions (SEDs) that depends more strongly 
on to their stellar mass than their star formation activity, and also 
is much less affected by dust extinction effects. 

Therefore, we decided to combine these three surveys to 
compile a sample from which to extract our sources. We deter- 
mined the objects in common to the UDS DR3+ K-band and 
the SUBARU catalog from Simpson et al. (2010), and cross- 
correlated this list of objects with the SWIRE data in the XMM- 
LSS region using the RA and DEC positions to within an error 
of 1 arcsec. The final area covered is 0.77 deg 2 (see Fig. [TJ. In 
the borders of the Fig.Q] the deep is lower but its means that we 
can lose some Z-ERO/OG there (1 or 2). Star-galaxy separation 
was made with the Bz-K diagram: B—i' and B-K s colors and the 
FWHM of objects following the criterion by Daddi et al. (2004). 
In what follows, colors are measured in a 3 arcsec aperture, while 
magnitudes are total magnitudes from SExtractor MAG AUTO 
parameter (Bertin & Arnouts 1996). We adopted for each object 
the same aperture among all passbands to ensure that the colors 
were measured correctly (Miyazaki et al. 2003; Furusawa et al. 
2008). 

The next step was to calculate the photometric redshifts 
(zphot) using the usual codes from the literature. In our case, 
we have made trials with several ones, and in the end we have 
selected the LePHARE program (Ilbert et al. 2006) based on a 
SEDs fitting method, because it gave the best results in compar- 
ison with available XMM-LSS objects spectroscopic redshifts. 
These photometric redshifts were computed with a large set of 
templates (Polletta et al. 2007), covering a broad domain in pa- 
rameter space with a Calzetti & Heckman (1999) extinction law. 

We have also compared our catalog with other zphot stud- 
ies in the literature in the UDS field. Williams et al. (2009) have 
a large sample of objects with zphot calculated with the EAZY 
code (Brammer et al. 2008) and Rowan-Robinson et al. 2009 had 
used its own code, IMPZ code (Rowan-Robinson et al. 2002). 
All these codes have reasonable sources of error and this behav- 
ior is more evident at higher z. Therefore, we have estimated a 
standard deviation (<x ; ) of four catalogs (in some sources only 
three redshifts are available): the Williams et al. catalog, the 
Rowan-Robinson et al. catalog and our results using apertures 
of 2" and 3" (our standard value in this paper). 

Among these sources, we consider only those for which 
z > 2.5 that satisfy the selection criteria of Z-EROs/OGs given 
in section|2]by the Eq. ((5]). There are not many available sources 
that comply with all of our constraints at z > 2.5. We find only 
45 galaxies, which decreases to 20 objects (see Fig.[T]and Table 
[1} when we select only those for which cr z /(l + z) < 0.1. That 
is, we reject the galaxies whose photometric redshift appears to 
depend on either the aperture size or the algorithm of photomet- 
ric redshift determination. This is the sample that we use for our 
analysis to ensure that we are conservative in the redshift range 
2.5 < z < 3.8. 

As noted by L10, any compilation of EROs/OGs is not a 
homogeneous sample of galaxies with similar intrinsic charac- 
teristics at all redshifts. The term 'Z-ERO' reflects the observed 
characteristic colors of a galaxy, not its intrinsic properties, and 




Fig. 1. Positions of the 63,550 galaxies of the merged catalogs 
SUBARU, UKIDSS, and Spitzer. The light circles correspond 
to the 20 selected galaxies complying with the criteria for Z- 
EROs/OGs given in section|2] 

Table 1. General characteristics of the final 20 Z-EROs/OGs 
(objects with asterisks represent the galaxies with <x,/(l + z) < 
0.05.) 



UDS Id 


RA (deg.) 


DEC (deg.) 


z 


(B-V),.,„ 


tabs (Gyr) 


06165 


34.43191 


-5.48688 


2.511 


0.786 


2.553 


79102 


34.14750 


-4.77846 


2.535 


0.920 


2.355 


25371 


34.73086 


-5.29274 


2.536 


0.834 


2.238 


01559 


34.40150 


-5.53476 


2.540 


0.655 


2.511 


40812* 


34.17666 


-5.13515 


2.565 


0.771 


2.716 


08303* 


34.29894 


-5.46366 


2.581 


0.770 


2.536 


24756 


34.73077 


-5.29849 


2.668 


0.797 


2.665 


65143 


34.40916 


-4.91072 


2.688 


0.904 


2.173 


78891 


34.59163 


-4.78038 


2.705 


0.583 


2.360 


51547 


34.90864 


-5.03568 


2.729 


0.689 


2.684 


39991* 


34.73174 


-5.14295 


2.754 


0.888 


2.574 


11488 


34.64674 


-5.43199 


2.813 


0.868 


2.712 


85772 


34.38530 


-4.71782 


2.852 


0.583 


2.487 


28153 


34.28822 


-5.26259 


2.889 


0.395 


2.043 


28841* 


34.14779 


-5.25500 


2.896 


0.798 


2.716 


44569 


34.42194 


-5.10082 


2.948 


0.544 


1.758 


10585 


34.79973 


-5.44024 


3.033 


1.060 


2.394 


26888 


34.24049 


-5.27663 


3.112 


0.658 


2.430 


14598* 


34.68966 


-5.40086 


3.287 


0.680 


2.744 


93574* 


34.28374 


-4.63816 


3.744 


0.705 


2.309 



the term 'OG' denotes that the object is an old early-type galaxy 
rather than a starburst. In any case, the range of stellar masses, 
M/L ratios, etc. may vary significantly among the 20 galaxies, 
and there may be selection effects. Nevertheless our analysis dis- 
tinguishes the dependence on redshift from that on luminosity. 



4. The rest-frame color (B - V). 

We use AB apparent magnitudes that are corrected for Galactic 
extinction (assuming the average value for the field of Ay=0.07), 
for different wavelengths, niABiM), = 1, ■•■,^V/), for 5 < Nf < 
12, with the corresponding error bars (in our case in the optical, 
near-infrared and mid-infrared, between B and [8.0]; we note 
that the [24.0] filter is not used here). As we have already men- 
tioned, we consider only data points with a flux S/N above 3. 
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Table 2. The properties of the final 20 objects with cr z /(l + z) < 0.1 subdivided into three bins of N galaxies. 



<f obs ) (Gyr) 


N 


{(B-V) Test ) 


(L v ) (1O 1U L0) 


<M„> (10 1U MO) 


(t g a) (Gyr) 


(fform) (Gyr) 


2.18 
2.50 
2.70 


1 
7 
6 


0.698 ± 0.200 
0.772 ±0.163 
0.768 ± 0.072 


14.0 ±4.9 
11.7 ±6.8 
8.0 ±6.5 


9.3 
11.3 
7.5 


i.3o 

1.82 + lf, 

i.8o il 


0.88 + V 2 i 
0.68 + ™ 
90 

U ' W -0.83 



Table 3. Properties of the six objects with crj{\ + z) < 0.05 (only one ? b s bin). 



<? obs > (Gyr) 


N <(B - V0 rest > 


(Lv) (1O 1U L0) 


<M„> (10 1U MO) 


<W ( G F) 


(ffo,m> (Gyr) 


2.33 


6 0.655 ±0.188 


16.5 ±5.5 


10.6 


1 20 +1 - 41 

i - zu -0.48 


i i q +0.48 

_ lj4 ] 



Applying the method of L10 (§3), we are able to determine the 
rest-frame color (B-V) of our selected galaxies. 

In Fig. |2 we display the average colors as a function of the 
age of the Universe when the galaxy is observed 

i r z -l l 

tobsXZ) = — dx- (7) 

Ho J°° i + x yjajl + x y + n A 

In this figure, we have also added the galaxies at z < 2.5 
from LI 0. The average gradient in color of the combined sample 
(including all galaxies in Fig. [3] 0.8 < z < 3.8) is 0.0877 + 
0.0093 Gyr 1 . For the galaxies at z > 2.5 analyzed in this paper, 
the gradient is 0.20 + 0.13 Gyr -1 . As noted by L10, the gap in 
the galaxy data points in the lower right corner of Fig. [3] is at 
least partially an artifact of the sample intrinsic color bias as a 
function of redshift. The step at z ~ 2.5 may be produced by 
the change in the color criterion to select EROs (z < 2.5) or Z- 
EROS (z > 2.5), bearing in mind that our criteria for selecting Z- 
EROs emulates the colors of EROs at z = 1 .5 at higher redshifts 
(see ©, hence we expect to find that the average color of Z- 
EROs is similar to those of EROS at z = 1.5. Therefore, one 
must interpret Fig. [3] as a plot of the intrinsic color versus (vs.) 
the redshift of the objects selected as (Z)-EROs/OGs, not as a 
general characterization of early-type galaxies. 

If we perform a bi-linear fit of the colors as a function of two 
independent variables t bs. and Merest, we get 

(B - V) resl =m+ a 2 [W(Gyr) - 2] + a 3 (M v<rest + 23), (8) 

with a x = 0.69 + 0.06, a 2 = -0.02 + 0.14, a 3 = 0.14 + 0.06 for 
the 20 galaxies in this paper. Given the large error bars, we are 
unable to discern any gradient of colors with respect to age and 
metallicity. L10 did measure significantly non-zero values of a 2 
and A3 at z < 2.5. 

5. Age estimation for early-type galaxies 

To estimate the average age and stellar mass of our galaxies, we 
apply the method described in L10 (§4). We assume a single stel- 
lar population for each galaxy; when there is instead a mixture 
of several populations, our ages would reflect an average of the 
age for the different populations. This method uses the synthesis 
model of Vazdekis et al. (1996), and breaks the age-metallicity 
degeneracy using the mass-metallicity correlation. We note that 
the role of thermally pulsing asymptotic giant branch (TP-AGB) 
stars is included in this synthesis model, and although some un- 
certainties are caused by this population (Maraston 2005), these 



are negligible for A < 6000 A (Bruzual 2007); these uncertain- 
ties are, however relevant for redder colors. Moreover, the ef- 
fect of TP-AGB stars is important for ages younger than 1 .5 Gyr 
and, as we later show, our galaxies are probably at this limit. 
For galaxies younger than 3 Gyr, the errors in metallicity are 
insignificant (L10), hence any evolution in the mass-metallicity 
relationship from high to low redshift would not significantly af- 
fect our results. 

Although the method cannot be applied to each single galaxy 
separately, it can be used statistically, by binning the galaxies 
within a small range of redshift and calculating the average age 
and its r.m.s. for the Z-EROs/OGs in each bin. In our case, we di- 
vided our sample into three bins, each of which has an associated 
average intrinsic color (see Fig. [3] and Table 0. The age of the 
early-type galaxies is plotted in Fig. [4] together with the ages of 
the galaxies with z < 2.5 from L10. The bin with 0.8 < z < 0.9 
is not plotted in Fig.|4]because it has a very large error bar (L10) 
and it does not provide any significant information. The vertical 
bars include the errors derived from the uncertainties in the in- 
trinsic colors plus the uncertainties in the synthesis model used 
in L10. We are neglect the systematic errors in the photomet- 
ric redshifts; were they not negligible, we would have additional 
systematic errors in the calculated ages. Figure [4] represents the 
average age of the given sample among the (Z-)ERO/OGs, with 
all the selection effects associated with each redshift. 

In Fig. [5] we plot the stellar masses in each bin derived with 
the same code. The errors in the stellar mass depend on both 
the data and variations in the initial mass function (IMF) in the 
Vazdekis et al. (1996) model, hence we can expect uncertainties 
of a factor two or even larger. The stellar masses of the galaxies 
at z > 2.5 are higher than 5 x 10 10 M Q . We note that at very 
high redshift, the stellar masses are much increased, as expected 
because of a selection effect that prevents us from detecting low 
mass galaxies at those distances. To differentiate between the 
effects of the dependence on mass and evolution, we apply a bi- 
linear fit to the average color, weighted with the square inverse 
of the relative errors in the average age, including the z < 2.5 
bins of L10. We get 

(B - V) resl =h+ b 2 (t obs , ~2) + b 3 log 10 (M.) , (9) 

fgai, = ci + c 2 (f ob s. - 2) + c 3 log 10 (M„) , (10) 

where both f ga i. and f bs. are in units of Gyr, M, is in units of 
10 11 M , and°/>, = 0.63 + 0.04, b 2 = 0.096 + 0.008, b 3 = 0.21 + 
0.06, ci = -0.18 + 0.96, c 2 = 1.90 + 0.29, and c 3 = 4.76 + 1.77. 
The average epoch of star formation (the first stars might be form 
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Fig. 2. Selection criteria used in this work to find the Z- 
EROs/OGs. a) Galaxies with z between 2.5 and 3.5. b) Galaxies 
with z between 3.5 and 4.5. The red open circles represent the 
final 45 candidates to Z-EROs/OGs following the Equation ((5]). 
The filled squares represent the final 20 selected galaxies with 
the zphot criterion. 



Fig. 3. (B — V) res t (Vega calibration) colors observed at different 
ages f b. s . (the age of the Universe corresponding to the given 
redshift). The XMM-LSS field are the data of this paper (the six 
filled triangles correspond to galaxies with < 0.05), ECDFS, 
and "Miyazaki et al." galaxies are taken from L10. 
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earlier) is ff 0rm . = f b s . - f ga i., as shown in Fig. [6] Separating the 
evolution from the mass dependence, 

fform. = di + d 2 (t obs . -2)+d } log 10 (M„,), (11) 

wheredi = 2.18+0.96,</ 2 = -0.91 +0.29, and d 3 = -4.76 + 1.77. 
For galaxies 0.8 < z < 2.5 L10 had found that d x = 1.94 ± 0.51, 
d 2 = -0.46 ± 0.32, and d 3 = -0.81 ± 0.98, hence with the ex- 
tension of the redshift to 3.8 in this paper, we have got with 
respect to L10 measure values of d 2 and do, that differ signifi- 
cantly from zero. In L10, the stellar masses at 0.8 < z < 2.5 
were confined to a narrow range, although for the wider red- 
shift range 0.8 < z < 3.8 the range of masses is similarly 
much larger (see Fig. 0, owing to Malmquist bias or other se- 
lection effects, permitting a more reliable determination of d$. 
The meaning of our numbers d 2 and dj, is that, the lower the 
redshift (for a fixed mass), the lower the age of galaxy forma- 
tion, and the higher the stellar mass (for a fixed redshift) the 
lower the age of galaxy formation. On average, for the whole 



Fig. 4. Average age (t ga i) of the stellar populations of (Z-)EROs. 
The dashed line stands for the limiting maximum t ga i - f Q b s . 
within the standard cosmology. The bin with 0.8 < z < 0.9 is not 
plotted because it has a very large error bar (L10) and it does not 
provide any significant information. 

range 0.8 < z < 3.8, we find that (fr^f = -0.91+0.29), 
( I to ™m ~ -4-76 + 1.77) Gyr, which is significant at the 3.1-cr 
and 2.1 -cr levels respectively. 

6. Discussion 

We have established a new color criteria for identifying objects 
of the same intrinsic colors (or even redder) than the EROs/OGs 
at redshift 0.8 < z < 2.5 but at higher redshift (Z-EROs/OGs). 
This is useful for searching for the oldest galaxies at an epoch 
during which the Universe was younger than ~ 3 Gyr. 
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Fig. 5. Average stellar mass of (Z-)EROs which are passively 
evolving early-type galaxies. 
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Fig. 6. Average age of the Universe at which the stellar popula- 
tions of the galaxies observed at age f b s . have been formed. 



By adding 20 galaxies at redshifts 2.5 < z < 3.8 to the cat- 
alog of L10, we have been able to see that we still derive for- 
mation ages that are within a narrow range (see Fig. [6), of on 
average 2.0 + 0.3 Gyr (zf orm . « 3 — 4) for the galaxies observed 
at 0.8 < z < 2.5 (L10), and 0.9+°* Gyr [simple weighted av- 
erage of the three bins in Fig. [6), and (zf 0lm . > 4.7, lcr) for the 
galaxies observed at 2.5 < z < 3.8. Clearly, the new galaxies 
observed at higher redshift were formed (around 1 Gyr) earlier, 
mainly because of the higher stellar mass. This is the most as- 
tonishing result derived in this paper: that very massive galaxies 
were formed at redshifts > 4.7. 

If we wish to be conservative, we could analyze only the 
galaxies with crz/(l+z) < 0.05. In this case, we would have only 
six galaxies and only one bin in our calculations (see Table [3] 
and the black triangles in Fig. [3). There is a slightly bluer aver- 
age color and higher average luminosity, but still the galaxies are 
older than 1 Gyr on average, and compatible with the results of 
the 20 galaxies within the error bars. Therefore, our results are 



quite insensitive to the level of robustness of the redshift deter- 
mination. 

This agrees with the results mentioned in the introduction 
that very massive evolved galaxies detected at redshifts 1.5-6 
were formed in the very early universe (Glazebrook et al. 2004, 
Daddi et al. 2005: Chen & Marzke 2004; Rodighiero et al. 2007; 
Wiklind et al. 2008). Clearly, any hierarchical A-CDM model 
that concludes that very massive galaxies formed after the for- 
mation of low-mass galaxies similarly disagrees with these re- 
sults. 

Acknowledgments: We thank A. Hempel (IAC) for help- 
ful suggestions and comments on a draft of this paper. 
We thank to A. Vazdekis (IAC) for clarifications about his 
synthesis model (Vazdekis et al. 1996). Thanks are also 
given to the anonymous referee for helpful comments. This 
work is based in part on data obtained as part of the 
UKIRT Infrared Deep Sky Survey; on observations with 
the Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology un- 
der a contract with NASA; and on data collected at Subaru 
Telescope, which is operated by the National Astronomical 
Observatory of Japan. NC is member of the Consolider- 
Ingenio 2010 Program grant CSD2006-00070: First Science 
with the GTC (http://www.iac.es/consolider-ingenio-gtc)" and 
the project EMIR AYA2009-06972 both of Spanish MICINN; 
and member of the SPIRE consortium led by Cardiff University. 



References 

Brammer, G. B., van Dokkum, P. G., & Coppi, P., 2008, ApJ, 686, 1503 
Bruzual, A. G, 2007, In: Stellar Populations as Building Blocks of Galaxies, A. 

Vazdekis, R. Peletier, Eds., Cambridge University Press, Cambridge, p. 125 
Calzetti, D., & Heckman, T. M. 1999, ApJ, 519, 27 

Casali, M., Adamson, A., Alves de Oliveira, C. et al. 2007, A&A, 467, 777 
Cassata, P., Cimatti, A., Kurk, J., et al. 2008, A&A 483, L39 
Chen, H.-W., & Marzke, R. O. 2004, ApJ 615, 603 
Cimatti, A., Cassata, P., Pozzetti, L., et al., 2008, A&A, 482, 21 
Daddi, E., Renzini, A., Pirzkal, N., et al. 2005, ApJ 626, 680 
Daddi, E., Cimatti, A., Renzini, A., et al. 2004, ApJ 617, 746 
Daddi, E., Cimatti, A., Pozzetti, L., et al. 2000, A&A, 361, 535 
De Lucia, G., Springel, V., White, S. D. M, Croton, D., & Kauffmann, G. 2006, 
MNRAS 366, 499 

Dickinson, M., Papovich, C, Ferguson, H. C. & Budavari, T. 2003, ApJ, 587, 25 

Fang, G.-W., Kong, X., & Wang, M. 2009, Res. Astron. Astrophys. 9, 59 

Fazio G. G., Hora, J. L., Allen, L. E. et al., 2004, ApJS, 154, 10 

Fontana, A., Pozzetti, L., Donnarumma, I. et al., 2004, A&A.424, 23 

Fontana, A., Santini. P., Grazian, A., et al. 2009, A&A, 501, 15 

Furusawa, H., Kosugi, G., Akiyama, M., Takata, T., Sekiguchi, K. & Furusawa, 

J., 2008, ASPC, 399, 131 
Georgakakis, A., Afonso, J., Hopkins, A. M., et al., 2005, ApJ, 620, 584 
Glazebrook, K., Abraham, R., McCarthy, P., et al., 2004, Nature, 430, 181 
Hempel, A., Cristobal-Hornillos, D., Prieto, M., et al, 2011, MNRAS, 414, 2246 
Hewett, P. C, Warren, S. J., Leggett, S. K. & Hodgkin, S. T., 2006, MNRAS, 

367, 454 

Ilbert, O., Arnouts, S., McCracken, H. J. et al., 2006, A&A, 457, 841 

Kriek, M., van der Wei, A., van Dokkum, P. G., Franx, M., & Illingworth, G. D. 

2008, ApJ, 682, 896 
Kriek, M., van Dokkum, P. G, Labbe, I., Franx, M., Illingworth, G. D., 

Marchesini, D., & Quadri, R. F. 2009, ApJ 700, 221 
Labbe, I., Franx, M., Rudnick, G., et al. 2007, ApJ 665, 944 
Lawrence, A., Warren, S. J., Almaini, O., et al., 2007.MNRAS, 379, 1599 
Lonsdale, C. J., Polletta, M., Surace, J.,et al. 2004, ApJS, 154, 54 
Lopez-Corredoira, M., 2010, AJ, 139, 540 (L10). 
Maraston, C, 2005, MNRAS, 362,799 
McCarthy, P., 2004, ARA&A, 42, 477 

Miyazaki, M., Shimasaku, K. & Kodama, T., et al. 2003, PASJ 55, 1079 
Miyazaki, S., Komiyama, Y., Sekiguchi, M., et al. , 2002, PASJ, 54, 833 
Polletta, M., Tajer, M., Maraschi, L. et al. 2007, ApJ, 663, 81 
Pozzetti, L. & Mannucci, F., 2000, MNRAS, 317, L17 
Rieke G. H, Young, E. T., Engelbracht, C. W., et al., 2004, ApJS, 154, 25 



Castro-Rodnguez & Lopez-Corredoira: Ages of red galaxies at high redshift 



7 



Rodighiero, G., Cimatti, A., Franceschini, A., Brusa, M., Fritz, J., & Bolzonella, 

M. 2007, A&A470, 21 
Roche, N. D., Almaini, O., Dunlop, J., et al., 2002, MNRAS, 337, 1282 
Rowan-Robinson, M., et al. 2008, MNRAS, 386, 697 
Rowan-Robinson, M., 2003, MNRAS, 345, 819 

Rudnick, G., Labbe, I., Forster Schreiber, N. M., et al. 2006, ApJ 650, 624 
Simpson, C, Almaini, O., Cirasuolo, M., et al., 2006, MNRAS, 373, L21 
Simpson, C, et al. 201 1, in prep. 

Somerville, R. S. 2005, in Multiwavelength Mapping of Galaxy Formation and 
Evolution, ed. R. Bender. & A. Renzini (Heidelberg: Springer,) 131 

Toft, S., van Dokkum, P., Franx, M., Thompson, R. I., Illingworth, G. D., 
Bouwens, R. J., & Kriek, M. 2005, ApJ 624, L9 

Van Dokkum, P. G., Quadri, R., Marchesini, D. et al. 2006, ApJ, 638L, 59 

Vazdekis, A., Casuso, E., peletier, R. F, & Beckman, J. E., 1996, ApJS, 106, 307 

Williams, R., et al. 2009, ApJ, 691, 1879 

White, S. D. M. & Rees, M. J., 1978, MNRAS, 183, 341 

Wiklind, T., Dickinson, M., Ferguson, H. C, Giavalisco, M., Mobasher, B., 

Grogin, N. A., & Panagia, N. 2008, ApJ, 676, 781 
Yan, Lin, Choi, Philip I., Fadda, D., Marleau, F. R. et al., 2004, ApJS, 154, 75 



